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ABSTRACT 

We have performed mid-infrared imaging of Barnard’s Star, one of the nearest stars to the 
Sun, using CanariCam on the 10.4 m Gran Telescopio Canarias. We aim to investigate an 
area within 1-10 arcsec separations, which for the 1.83 pc distance of the star translates to 
projected orbital separations of 1.8-18 au (P > 12 yr), which have not been explored yet with 
astrometry or radial velocity programs. It is therefore an opportunity to enter the domain of 
distances where most giant planets are expected to form. We performed deep imaging in the 
(V-band window (Si-2 filter, 8.7/jm) reaching a 3o detection limit of 0.85 ±0.18 mJy and 
angular resolution of 0.24 arcsec, close to the diffraction limit of the telescope at this wave¬ 
length. A total of 80 min on-source integration time data were collected and combined for the 
deepest image. We achieved a dynamical range of 8.0 ± 0.1 mag in the 8.7 /./tn band, at angu¬ 
lar separations from ~2 to 10 arcsec and of ~6-8 mag at 1-2 arcsec. No additional sources 
were found. Our detectability limits provide further constraints to the presence of substellar 
companions of the Barnard’s Star. According to solar metallicity evolutionary models, we can 
exclude companions of masses larger than 15 Mj up (T e g > 400 K), ages of a few Gyr, and lo¬ 
cated in ~3.6-18 au orbits with a 3O confidence level. This minimum mass is approximately 
5 Mj up smaller than any previous imaging survey that explored the surroundings of Barnard’s 
Star could restrict. 

Key words: brown dwarfs - stars: imaging - stars: individual: Barnard’s Star - solar neigh¬ 
bourhood - infrared: planetary systems. 


1 INTRODUCTION 

Optical and near-infrared high-contrast observations are of ma¬ 
jor importance in the study of substellar objects. Direct imaging 
searches for brown dwarfs and planets around stars explore a range 
of physical separations complementary to that of radial velocity 
or transit methods and provide key information on their forma¬ 
tion processes. Benefit of detecting the direct light of planets is 
that it enables their extensive characterization. We can determine 
their ages, masses, radii, effective temperature and, in particular. 


through spectroscopy we can have insight into their complex at¬ 
mospheres (Marley et al,|2012| |Currie et al.|20 13 ( i . Yet, reaching 
a sufficiently high-luminosity contrast at arcsecond level angular 
separations from the target star is still very challenging. As a result, 
in comparison with the transit method or precision radial velocity 
measurements, direct imaging brought only a few discoveries so 


far, e.g., jChauvin et al. (2004). Bejar et al. (2008 

i, Marois et al. 

(2008 

2010 ), Lagrange et al. 

(2010 

). Rameau et al. 

(2013), Carson 

et al. ( 

2013 ) and Gauza et al. 

(2015 

I. Most of the planets found by 


imaging are massive (M >5—10 Mj up ), young (x < 500 Myr), and 
are located at relatively large projected physical separations (~ 10— 
30 to 1000 au) around their host stars. Thus, each single discovery 
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Table 1. Observing log of GJ 699 with CanariCam at the GTC 


OB 

Observation 

date 

Start time 
(UTC) 

Save sets per 
nod beam 

Nod 

cycles 

Saveset 
time (s) 

On-source 
time (s) 

Instrument 

PA(°) 

PWV 

(mm) 

Readout 

mode 

1 

2012 July 29 

23:13:18.7 

8 

10 

6 

3x403.7 

0 

8.6-9.3 

S1RLCR 

2 

2012 July 30 

00:20:01.6 

8 

10 

6 

3x403.7 

90 

8.6-9.3 

S1RLCR 

3 

2013 June 09 

03:40:42.2 

29 

12 

1.5 

3x431.9 

0 

6.7 

S1R3 

4 

2013 June 10 

04:18:36.2 

29 

10 

1.5 

3x359.9 

300 

6.7 

S1R3 


extending the known population of imaged planets has still an im¬ 
portant impact on the field. 

In the formation process, gravitational collapse energy is re¬ 
leased and heats the interior of a planet but, as there is no internal 
source of energy, planet cools and fades down with age. Giant plan¬ 
ets at ages younger than 500 Myr are more easily detectable since at 
that stage their self-luminosity is still significant. However, the vast 
majority of stars from the solar vicinity (d ~2-10pc) are relatively 
old, with ages similar to the age of the Sun. Any potential planet 
from a nearby system, unless orbiting at close separation from its 
host star, where it is strongly irradiated, is expected to have cooled 
down to effective temperature below 500-600 K. Its spectral energy 
distribution would peak in the mid-IR range. Moreover, a solar-type 
star will be relatively faint at these wavelengths, so that the con¬ 
trast necessary to detect a planet or a brown dwarf companion will 
be lower. Currently, direct imaging of planets is feasible with the 
use of Hubble Space Telescope (HST) or the largest ground-based 
facilities equipped with adaptive optics systems, e.g. Gemini/NICI, 
Keck/NIRC2, VLT/NaCo, Subaru/HiCIAO, operating in the opti¬ 
cal or near-IR regime. Searches require proper observing and data 
processing techniques which attempt to remove diffracted light like 
Lyot coronograph y (|Sivaramakrishnan et al.|2001) > or nulling inter¬ 
ferometry ( |Hinz et al. ||1998|~ Serabyn 2000) and to suppress the 
speckle background (Angular Differential Imaging; |Marois et al.| 
2006; Lafreniere et al. 20071. In this work, we explore the potential 
of ground-based imaging at mid-IR wavelengths to directly detect 
planetary-mass objects. 

Barnard's Star (GJ699, V25000ph) is the fourth-closest in¬ 
dividual star and the second-closest system to the Sun currently 
known, after the triple aCentauri stellar system. It is also the 
closest star in the Northern hemisphere. It was classified as an 
M4.0V-type red dwarf l |Reidetal.|1995) , located at 1,824±0.005 pc 
l |van Le euwen 2007) and moving at the highest known proper mo¬ 
tion p = 10.37 arcsec per year. A slow rotation with a period of 
about 130 d ((Benedict et al. 1998), low magnetic activity (Htinsch 
|et al.|[l999) and other age indicators suggest that it is probably 
older than the Sun. Since its large proper motion was measured by 
|Bamard| l [l916) , the star received much attention from astronomers. 
Its bolometric luminosity and effective temperature were found to 
be (3.46±0.17) x 10~ 3 L s and 3134± 102 K (Dawson & De 


|Robertis|2004). From mass-luminosity relations for very low mass 
stars, |Delfosse et al,| (2000 ) determined the mass of GJ 699 to be 
0.158±0.013 M 0 . Particularly, extensive studies have been carried 
out refining the limits of possible planetary-mass companions (e.g., 
(Benedict et al.|1999l|Dieterich et al.|2012[|Choi et al.|2013) . Here, 
we present the results of mid-IR imaging of this star using Ca¬ 
nariCam at the GTC. The following sections of the paper describe 
the performed observations, reduction and analysis of the collected 
data, determination of the sensitivity limits, its translation to the 
physical parameters (mass, separations) of detectable companions 
and comparison with the previous surveys. 



Figure 1. Final image of the Barnard’s Star at 8.7 pm obtained with Canari¬ 
Cam at the GTC. The used on-source time is 80 min. Field of view shown 
in the left-hand image is 13 x 13 arcsec. Right-hand image is the enlarged 
5x5 arcsec central part with a different contrast scale, marked with a grey 
square on the larger field image. North is up, east to the left. 


2 CANARICAM OBSERVATIONS 

Observations of the Barnard’s Star were carried out in queue mode, 
during the nights of 2012 July 29 and 2013 June 9 and 10 (UT). 
We used the mid-infrared camera CanariCam (Telesco et al.|2 008) 
operating at the Nasmyth-A focal station of the 10.4 m Gran Tele- 
scopio Canarias (GTC) at the Roque de los Muchachos Obser¬ 
vatory on the island of La Palma. CanariCam was designed to 
reach the diffraction limit of the GTC at mid-IR wavelengths (7.5- 
25 /rm). The instrument uses a Raytheon 320x240 Si:As detec¬ 
tor with a pixel scale of 79.8 ±0.2 mas, which covers a field of 
view of 25.6 x 19.2 arcsec on the sky. We imaged our target in the 
10 micron window, using a medium-band silicate filter centred at 
X — 8.7 pm (AX =1.1 /rm). The choice of this particular bandpass 
was a compromise between the instrument performance, in partic¬ 
ular the filters transmissivity, and the sky background contribution 
significantly higher at the N broad-band and other narrow-band fil¬ 
ters than at the Si-2 filter. Si-2 is also favoured by a better spa¬ 
tial resolution, since the diffraction disc is larger at the available 
narrow-band filters at longer wavelengths. 

Observations were performed with the standard chopping and 
nodding technique used in the mid-IR to remove the sky emission 
and radiative offset. Chopping consists of switching the telescope 
secondary mirror at a typical frequency of a few (2-5) Hz be¬ 
tween the position of the source (on-source) and the nearby sky 
(off-source). This rapid movement of the secondary mirror allows 
subtraction of the sky background emission that is varying in time 
at frequencies below the chop frequency. Movement of the sec¬ 
ondary mirror changes the optical configuration of the telescope, 
resulting in two different emission patterns seen by the camera and 
producing a spurious signal termed the radiative offset in the chop- 
differenced images. To remove the radiative offset, the telescope is 
moved between two nod positions to swap over on- and off-source 
positions. We used a ABBA nodding sequence and ‘on-chip’ chop- 
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Figure 2. Left: image mosaic of the Barnard’s Star, processed and stacked to encompass the full area available in the CanariCam data. Right: WISE W3 image 
of the Barnard’s Star, after subtraction of the PSF of the target. The yellow circle on both images marks an area within 14 arcsec radius around the centroid 
position of the star. The CanariCam mosaic image fully covers the area up to this separation and extends up to 20 arcsec but with incomplete coverage and 
lower sensitivity. The shown field of view of the WISE image cut-off is 70x70 arcsec. No additional source is detected within this separation. Residuals of the 
star after PSF subtraction extend to roughly 14 arcsec. North is up and east is to the left. 


ping and nodding, with a chop-throw and nod offset of 8 arcsec, a 
chopping frequency of 1.93-2.01 Hz and a nod settle time of about 
45 s. On chip method is recommended whenever the scientific tar¬ 
get is point-like, since both on-source and off-source chop positions 
contain the signal of the target inside the detector field of view and 
can be aligned and combined. On-source integration time in each of 
the four observing blocks (OBs) was 20 min, divided into three data 
files composed of a set of images (savesets ) at subsequent chopping 
and nodding positions. Individual frames of 26 and 19 ms expo¬ 
sures were co-added by CanariCam control software to savesets of 
1.6 and 6 s for S1R3 and S1R1-CR readout modes, respectively. A 
detailed observing log is given in Table [T] 

In total, we integrated the source for 80 min. Sky condi¬ 
tions during the observations were photometric. Precipitable water 
vapour (PWV) as measured by the Instituto de Astroftsica de Ca- 
narias real-time PWV monitor, was 8.6-9.3 and 6.7 mm for OB 1, 2 
and OB 3, 4, respectively. We observed with a set of three different 
position angles of the instrument on the sky, 0°, 90° and 300°, and 
with two chop position angles of the secondary mirror, 0° and 90°, 
to guarantee that the chop offset was done along the X-axis. By 
that, we wanted to avoid the loss of potential objects in the regions 
overlapping with the negative images of the star from the off-source 
chop positions or in parts obscured by the horizontal cross-talk fea¬ 
tures that appear for bright sources jOkamoto et al.|2003) . Also, we 
could discard eventual contaminants along the chop axis. 


3 DATA REDUCTION AND ANALYSIS 
3.1 Image processing 

The data were processed using standard routines within the IRAlQ 
environment. CanariCam images are stored in the standard multi- 

1 IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in 
Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 


extension fits files, with a structure of [320,240,2, M][N], where 
320 and 240 are the image pixel dimensions, 2 is the number of 
chop positions, M of savesets and N of nod positions. 

Off-source savesets, where the position of the secondary mir¬ 
ror is not aligned with the primary mirror, were subtracted from 
the corresponding on-source savesets, for respective nod beam po¬ 
sition. These chop-/sky-subtracted frames, where the star is located 
at the centre of the detector, were then aligned to correct for very 
small misalignments (of less than 3 pixels), and each pair corre¬ 
sponding to the A and B nod positions were combined, to subtract 
the radiative offset. The sky-subtracted frames where multiplied 
by —1 to recover the negative contributions of the star (off-source 
position of the secondary mirror). Because the negatives in the A 
and B nod positions do not overlap, being at opposite sides and 
at 8 arcsec of the on-source central location, they were radiative- 
offset subtracted before they were aligned. Residual detector levels 
constant along single columns or lines but varying across these re¬ 
mained in both the positive and negative chop- and nod-subtracted 
frames; these were background fitted (masking the target) and sub¬ 
tracted. The alignment itself was applied at once, to all (positive 
and negative) images of consecutive repetitions of an OB, relative 
to a same reference image and so that the target has its centroid 
located on an integer pixel position. Before aligning, the images 
were copied into larger ones to avoid the trimming of outer data re¬ 
gions. Then the frames were average-combined per repetition or al¬ 
together, using a shallow sigma upper and lower clipping to discard 
occasional short transients and extreme pixel values. Each combi¬ 
nation involved masking the negative counts of the target. OB 4 
had cross-talk features that we had removed. OB 2 and 4 were ac¬ 
quired with position angles differing from the North-up East-left 
orientation. For OB 2, we transposed the stack to North-up East- 
left, whereas for OB 4, we resampled the stack. In total, the data of 
OB 1, 2 and 3 were thus resampled only once, whereas those of OB 
4 twice. For the combination of the stacks of the different OBs, the 
stacks were flux-scaled according to their zero-point magnitude - 
as measured on the target - and weighted inversely proportional to 
the scaled variance of their background noise and the square of the 
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Figure 3. The Si-2 band (8.7 pm ] 3 a contrast curves, derived on the deepest 
obtained image from artificial source injection (blue line with circle marks) 
and from background noise computed in concentric annuli around the star 
(red line with cross marks). We reach ASi-2 ~ 8 mag at separations larger 
than 2 arcsec. The dashed line shows the sensitivity of combined single 
epoch images OB 3+4, which reaches a detection limit of about 0.25 mag 
less beyond 2 arcsec separation. The two bumps at about 0.3 and 0.8 arcsec 
are caused by the presence of first and second Airy rings. 

full-width half-maximum (FWHM) of the target. The central part 
of the final reduced image, with total on-source time of 80 min is 
displayed in Fig. IT] An image mosaic of the whole covered area is 
displayed in Fig. [2] This all-epoch mosaic extends to distance of 
about 20 arcsec (37 au) to the target at the four cardinal directions. 


3.2 Sensitivity and detection limits 

The star point-spread function (PSF) on the final obtained images 
has an FWHM of 3.01 pixels corresponding to 0.24 arcsec, for the 
image including all-epoch data, and 2.89 pixels (0.23 arcsec) for 
the image including single epoch data from OB 3 and 4. It is close 
to the theoretical FWHM of the diffraction-limited PSF, which for 
GTC is 0.19 arcsec at 8.7 pm. 

To determine the Barnard’s Star magnitude in the Si-2 filter 
we used the J,H,K photometry from the Two Micron All Sky Sur¬ 
vey (2MASS; |Skrutskie et al.|2006| and W 1, W2, W3,W4 from the 
Wide-field Infrared Survey Explorer ( WISE) All-Sky and A11WISE 
Source Catalogs ( [Wright et al.|2010) . We converted the 2MASS 
and WISE magnitudes into fluxes using the corresponding Vega 
zero points and interpolated the values to obtain average flux at 
8.7 pm through a least-squares fit to a power function. WISE W2 
measurement was not used in the fit because of saturation, affect¬ 
ing sources brighter than approximately 6.7 mag in this band. The 
calculated Si-2 brightness of Barnard’s Star is 4.12 ± 0.19 mag, 
using the Vega system zero point determined for this CanariCam 
filter. This value is very close to the VP3 magnitude of the star 
(W3 = 4.036 ±0.016 mag). For our deepest image, that is, the one 
combined from all OBs, the limiting magnitude is 11.92 ±0.25 mag 
(0.85 ±0.18 mjy), estimated using the ratio of the peak counts of 
the star to 3c background noise. 

In Fig. [3] we present the 3cr contrast curves for CanariCam im¬ 
age of Barnard’s Star. The two solid lines shown in the plot corre¬ 
spond to two different contrast calculation methods that have been 
used. In the first approach, we computed the background noise, a, 
as a function of radial separation from the target star, by measuring 



Figure 4. Same as Fig. 1. but with artificial sources inserted along the di¬ 
agonals in left-hand panel, and after PSF subtraction and the addition of 
artificial sources in right-hand panel. 


the standard deviation in a 1 pixel wide concentric annuli around 
the star. The 3a noise counts were converted to contrast (delta mag¬ 
nitudes between the primary star and the measured quantity, noise 
in this case) by relating to the peak pixel value of the star PSF. The 
contrast curve obtained using this method is plotted in Fig. [3] with 
a red solid line with cross marks. 

In the second method, we used artificial sources to estimate 
the detection limits and contrast as a function of radius. Artifi¬ 
cial objects were produced using the PSF model of the primary 
and inserted on the diagonal axes, with intervals of ~ 2 x FWHM 
(~ 0.6 arcsec) and decreasing steps in the inner parts of the image. 
We considered as a detection when the added object at a given mag¬ 
nitude and radial distance is detected through visual inspection in at 
least three of the four diagonal positions. A simulated source with 
its PSF retaining characteristic stellar shape and visible marginally 
above the background noise corresponds typically to a signal to 
noise of 3-5. To examine the inner area (within 1 arcsec) around 
the star, we performed target PSF subtraction on the images with 
inserted artificial sources by rotation or flipping along an axis, by 
previously measuring the centre and semi-major axis of the ellipse 
isophote closest to the level of the candidates (~10c) we expect 
to find in the PSF wings. Example images with simulated artificial 
objects are shown in Fig. [4] Results from this two methods, which 
were found consistent, yield that we reach a dynamical range in Si- 

2 of about 6 mag at ~ 1 arcsec separation from the target star and a 
maximum of 8 mag at separations >2 arcsec. The second method 
systematically yields detection limit as a function of projected sep¬ 
aration from the Barnard’s Star that is about 0.25 mag brighter than 
obtained with the first method. 

At the range of projected orbital separations explored in this 
search, with roughly one year period separating the two observa¬ 
tions, the orbital motion of a potential companion may not be neg¬ 
ligible. We estimate that a planetary mass companion at a 18 au 
circular orbit orientated face-on would move ~0.3 arcsec over one 
year baseline which is approximately the FWHM of our CanariCam 
images, and up to 0.8 arcsec at 3 au orbit (1.5 arcsec angular sep¬ 
aration). At an edge-on orientation, the shift could be higher than 
0.3 arcsec at separations closer than 5 au. Since for some range of 
parameter space the displacement would be higher than the spa¬ 
tial resolution of the images, we provide also the detection limit of 
single epoch data. The dashed line in Fig.|3]is the contrast curve de¬ 
termined on the image obtained using only first epoch images (OBs 

3 and 4, TablejTjl. The achieved detection limit in this case, with a 
total of 40 min on-source time, is about 0.25 mag lower than that 
achieved with the use of both epochs observations. 
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Mass [M Jup ] Mass [M Jup ] 

Figure 5. Theoretical absolute magnitudes versus mass of brown dwarfs and giant planets at Si-2 8.7 /jm (left) and N bands (right). For Si-2 we used the 
Ames-COND models at ages of 0.1, 0.5, 1, 5 and 10 Gyr and for A-band we used models by |Burrows et al. j(l997). Dashed horizontal lines mark our 30 
detection limit range from single epoch observations: Me,, 2 = 15.36±0.28 and from the deepest image obtained using both epochs: Msi _2 = 15.62±0.25 mag. 
Also several isotherms are plotted in the left-hand panel with grey dotted lines. According to the models, this implies that companions more massive than 
~ 15 Mj up or with T c ff > 400 K at solar ages would have been detected. 


4 CONSTRAINTS ON THE PRESENCE OF 
COMPANIONS 

4.1 Physical interpretation of detection limits 

In the following we relate the detection limits of our search to 
the physical properties of brown dwarf and planetary companions, 
specifically to their masses, effective temperatures and luminosi¬ 
ties. A complication inherent to objects in the substellar domain is 
the continuous cooling in the course of their evolution. It precludes 
the possibility to estimate the mass applying unique relations inde¬ 
pendent of age, such as the mass-luminosity relation for the main- 
sequence stars. In this case, we need to rely on theoretical models 
providing a grid of luminosities, temperatures, synthetic photom¬ 
etry as a function of masses and ages. In this work, we used the 
Ames-COND models dAllard et al.|200ll|Baraffe et al.|2003||ATl 
|lard et al.|2012[ l for solar metallicity and the models of giant plan¬ 
ets and brown dwarfs by |Burrows et al.| l | 1 997} . Both the COND 
and the Burrows et al. (1997) models apply to T e g- < 1300 K and 
extend down to 100 K. They include the formation of dust in the 
atmospheres of this objects, however they neglect its opacity, con¬ 
sidering that the dust grains settle below the photosphere. To com¬ 
pute the synthetic magnitudes for the Si-2 8.7 /jm band we used the 
Phoenix Star, Brown Dwarf & Planet Simulator available onlin^] 
We input the transmission file of the Si-2 filter and obtained the 
isochrones for a set of five different ages: 0.1, 0.5, 1, 5 and 10 Gyr. 

At the relatively old ages (> 1 Gyr), substellar companions in 
the planetary-mass regime have cooled down to temperatures below 
600 K, according to theoretical evolutionary models (Burrows et al. 
|1997| |Chabrier et al.|2 000 ). At these temperatures, brown dwarfs 
and giant planets are very faint in the near-IR and emit most of 
their flux in the mid-IR. As an example, while a 5 Mj up planet near- 
IR flux decrease from M] ~ 15 mag at lOMyr to Mj ~ 25 mag at 


- http://phoenix.ens-lyon.fr/simulator/index.faces 


4.5 Gyr, the mid-IR emission only changes from ~ 11.5 mag to 
Mn ~ 16 mag in the same interval of age. 

In the left-hand panel of Fig. [5] we plot the Si-2 absolute mag¬ 
nitudes versus masses produced by the models, in comparison with 
the 30 detection limit of our CanariCam observations. For ages of 1 
Gyr and younger, objects down to 5 Mj up would have been detected. 
At solar ages, we are sensitive to objects with masses higher than 
approximately 15 Mj up or with effective temperatures above 400 K 
and at projected separations from ~3 to 18 au. At smaller separa¬ 
tions, around 1 arcsec (~2 au), we could have detected objects with 
masses above ~20Mj up . 


Several indicators point that the Barnard’s Star belongs to the 
older population. The non-detection of lithium in the atmosphere 
imposes a minimum age of 20Myr (Zuckerman & Song 20041. 
The level of chromospheric activity suggest an age above 600- 
800 Myr ( jRauscher & Marcy|2006| l. It has a very low X-ray lumi¬ 
nosity (log L x = 26) indicating low level magnetic activity lJHunsch| 
|et al.|[]~999] |Vaiana et al.||198l) . |Eggen| j 1 996) estimated the age 
of GJ 699 at 10 Gyr based on the Ca II index. Also its high space 
velocities (Leggett 1992; Eggen 1996), lower than solar metallic¬ 
ity [Fe/H]=—0.39 ±0.17 Gizis|1997 ' and a probable long rotation 
period of approximately 130 days (Benedict et al. 19981 are all con¬ 
sistent with a relatively advanced age of 7-12 Gyr. At that age, our 
sensitivity limit enables us to detect companions more massive than 
approximately 20Mj up with r e g- > 450 K. 


For a comparison, we have checked the models by Burrows 
|et al.| ( |1997| right-hand panel of Fig. [5} that provide synthetic pho¬ 
tometry in the N band (10 /./m) at ages of 0.1, 0.5, 1 and 5 Gyr for 
objects with masses below 40Mj up . We find that the two models 
are in a fairly good agreement and give a similar mass and T e g con¬ 
straint using the same Si-2 sensitivity limit. We interpret that the 
slightly brighter magnitudes for the same masses yield by the N- 
band isochrones result from a wider bandpass and longer central 
wavelength of this filter, at which the flux of very cool objects is 
expected to be larger, as well as from differences between the two 
models. 
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4.2 Comparison with previous searches and WISE data 

First claim of planetary companions of Barnard’s Star were re¬ 
ported in the late 60s by |van de Kamp| ( |1969| l. Using astrometry, 
he detected perturbations in the proper motion of the star consis¬ 
tent with two planets comparable in mass with Jupiter. These plan¬ 
ets however were definitively ruled out by the radial velocity stud¬ 
ies. jCtoTeTar] (j20l3j used Doppler measurements obtained from 
Lick and Keck Observatories with a precision of 2 ms~* during 
an 8 yr monitoring to preclude planetary companions of masses 
above 2 Earth masses with periods below 10 days and above 10 
Earth masses with periods up to 2 yr, save for face-on orbits, (area 
labelled T in Fig. [6j. This result is supported by the probability 
of 94.6-98.2% that the orbital inclination of putative companions 
is larger than 11°-19°. Various ground- and space-based imag¬ 
ing surveys dedicated to detect ultracool companions targeted the 
Barnard’s Star along the past two decades. In this section, we fo¬ 
cus on the ones that reached the highest sensitivities up to date and 
compare them with our results. 

|Dieterich et al.| < |2012| > carried out a large, volume-limited 
search for substellar companions to stars (mainly around M dwarfs) 
in the solar neighbourhood within ~10 pc. The survey was per¬ 
formed with the HST NICMOS instrument, using four filters cen¬ 
tred at 1.10, 1.80, 2.07 and 2.22 /urn, for snapshot high-resolution 
imaging of 255 individual stars. Physical separations accessible 
in this search correspond to mean semimajor axes between 5 and 
70 au. The minimum detectable mass was estimated to be 42Mj up 
at 3 Gyr or 52Mj up at 5 Gyr. The masses were derived using the 
|Chabrier et al.|{2000| > models. For the Barnard’s Star however, as it 
is one of the nearest stars from the sample, closer separations were 
probed. For the survey limits on this particular star, we used the 
information from the table 2 in Dieterich et al. (2012), which gives 
the sensitivities achieved for each target at various separations. The 
range of masses and separations of companions ruled out by this 
survey is marked as area II in Fig. [6] 

|Oppenheimer et al~| < ]2()01 ft observed northern stars from the 8 
pc sample including GJ 699, in the optical (Gunn r and z filters) 
and near-infrared wavelengths (J,K filters). They used the Adap¬ 
tive Optics Coronograph instrument on the Palomar 1.5 m telescope 
for the optical imaging and the Cassegrain Infrared Camera on the 
Palomar 5 m Hale Telescope, for the near-IR. For about 80% of the 
surveyed stars, companions more massive than 40Mj up at 5 Gyr 
age would have been detected at separations between 40 and 120 
au. The mass detection limits were determined using the Burrows 
|et al.|{T997| evolutionary models. The area marked as III in Fig. [6 
shows the limits of the lowest mass detectable by the survey at the 
closest accessible separations, given in tables 11 and 12 in |Oppem | 
|heimer et al.|{200T} . 

|Schroeder et al.| {2000 1 searched for faint companions to 23 
stars within 13 pc of the Sun using the HST Wide Field Planetary 
Camera 2. Barnard’s Star was imaged through F615W, F8141V and 
F1042M filters using short and long exposures ranging from 1.6 to 
600 s, to map different separations from the target star. According 
to the limits shown in fig. 13 in the paper, their survey was sensi¬ 
tive to objects of about 40 A/j up at angular separations larger than 
1.5 arcsec, 30 Afj up at > 2.0 arcsec and 20 Mj up beyond 3.5 arcsec, 
for an age of 5 Gyr. The authors also employed the |Burrows et al.| 
{1997} models to derive the limits of the survey in terms of com¬ 
panion masses. These limits are shown in Fig. [6] as area IV, which 
extends up to ~30 au. 

Fig .[^{summarizes in a simplified manner the detectability lim¬ 
its of the described surveys and of the CanariCam imaging search in 



Figure 6. Constraints on the presence of substellar companions around the 
Barnard’s Star. Comparison of obtained CanariCam detection limits (thick 
line, shaded area) with results from previous searches and limits determined 
using WISE W3 image. Corresponding surveys marked in the diagram are 
described in Section 4.2. The CanariCam area with dashed line border cor¬ 
responds to the region at angular separations >14 arcsec, where the mosaic 
coverage is not complete. 

terms of companion mass lower limits and projected physical sep¬ 
arations. Additionally, we have included the detection limits of the 
search for companions we have carried using the WISE Wi band 
data at 12 fjm. We have retrieved Barnard’s Star’s image from the 
WISE All-Sky release using the IRSA Wise Image Service. To de¬ 
rive the sensitivity in the inner part, we have previously performed a 
PSF subtraction of the target using as a reference the PSF of a sim¬ 
ilar brightness star located close to the Barnard’s Star on the sky, 
after scaling their peak fluxes. We then obtained the 3ct limit of 
sensitivity of IV3 using the same method as for the CanariCam data 
and compared to the Ames-COND models computed for this band. 
We estimated a sensitivity of W 3 = 9.7-9.8 mag at 14—15 arcsec 
separations, 10.5 mag at 20 arcsec and a limit of 11.5 mag (~1.0 
mjy) at 40 arcsec and beyond. According to models, a (solar age) 
companion of at least 15 Mj up located at separation > 40 arcsec (> 
73 au) would have been detected. The obtained detection limit in 
terms of sensitivity and of objects masses and effective tempera¬ 
tures is in agreement with the general documentation of the mis¬ 
sion {Wright et al.|[20l0| ) and other studies that employed WISE 
data, e.g. jLuhman|{2014j >. This search for companions using WISE 
W3 data which extends to wider separations is complementary with 
our higher resolution CanariCam images and provides similar sen¬ 
sitivity. 

Radial velocity and astrometry techniques exclude planetary 
companions of less than 1 Mj up with orbital periods up to 2 yr 
{Choi et al.|20 1 3[ |Benedict et al.|199~9| >. This detection limits are 
still far beyond the capabilities of any other method, whereas at 
wider orbits direct imaging provides the strongest constraints on 
the presence of companions. With CanariCam, we could have de¬ 
tected companions more massive than 15 AJj up at projected separa¬ 
tions from ~3 to 18 au for an age of 5 Gyr. CanariCam observations 
have enabled us to set the strongest constraints on the presence of 
very low mass brown dwarfs in wide orbits (3.6-18 au) with a 99% 
confidence level. However, it should be noted that this statement is 
based on model evolution predictions. Were the models prove to be 
invalid for these least massive substellar objects, the obtained mass 
and temperature lower limits would have to be reviewed. In particu- 
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lar, the Ames-COND isochrones that were used, are available only 
for solar abundance. Given the slightly sub-solar metallicity of the 
Barnard’s Star, this may result in a difference of the derived mass. 


5 SUMMARY AND FINAL REMARKS 

We performed a deep, high spatial resolution imaging of the 
Barnard’s Star at mid-IR Si-2 8.7 ,um wavelength using Canari¬ 
Cam at the 10.4 m GTC telescope. No companion candidates 
were found on the obtained images. Our detectability limits pro¬ 
vide further constraints on the presence of substellar companions. 
With 80 min on-source integration time we achieved sensitivity 
of 0.85 ±0.18 mjy allowing to detect massive planets and brown 
dwarf companions down to 15 AJj up which corresponds to effec¬ 
tive temperatures above 400 K, assuming a solar age. Our search 
covers a field of l-10arcsec radius around the target star (~2-18 
au for Barnard’s Star distance), which means that we can probe the 
domain of distances where most giant planets are expected to form 
l |Chabrier et al.|2014| l. 

This work demonstrates that the modern ground-based 
mid-IR imaging instruments operating on 10-m class telescopes 
can reach angular resolutions and sensitivity limits as good and, in 
certain cases (e.g. nearby, relatively old stars) better than adaptive 
optics systems in the optical or near-IR or space telescopes. This 
technique presents a high potential to perform direct imaging 
studies of brown dwarfs and exoplanets. 
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